Attosecond streaking experiments on atoms: quantum theory versus simple model 
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A new theoretical approach to the description of the attosecond streaking measurements of atomic 
photoionization is presented. It is a fully quantum mechanical description based on numerical solving 
of the time-dependent Schrodinger equation which includes the atomic field as well as the fields of 
the XUV and IR pulses. Also a simple semiempirical description based on sudden approximation is 
suggested which agrees very well with the exact solution. 



One of the problems in the physics of ultrashort atomic 
phenomena is the characterization of the attosecond ex- 
treme ultraviolet (XUV) pulses which are used in pump- 
probe experiments studying the time evolution of atomic 
processes. The characterization includes measurements 
of the pulse duration, intensity, basic frequency, spectral 
distribution and possible chirp. In a typical experiment 
H H Q- the XUV pulses with a duration of a few 
hundred attoseconds are produced by filtering out the 
high energy part of the high harmonic spectrum which 
results from a nonlinear interaction of a few-cycle in- 
tense infrared (IR) laser pulse with a gas target. The 
laser pulse duration is typically 5-7 fs, its carrier wave 
length is ~ 800 nm (the photon energy El ~ 1.6 eV), 
thus the period of the laser pulse is larger than the XUV 
pulse duration. The intensity of the laser is typically 
3 x 10 13 — 10 14 W/cm 2 . At such intensity the electric 
field of the laser, Et, ~ (0.1 — 0.3) x 10 9 V/cm, is much 
smaller than the atomic field (1 a.u.= 5.14 xlO 9 V/cm) 
although it is strong enough to produce high harmonics. 
The electric field of the XUV pulse is still smaller, about 
7 x 10 7 V/cm 0, so that the linear response approxima- 
tion can be readily applied to the interaction of the XUV 
pulse with atoms. The typical carrier frequency of the 
XUV pulse corresponds to the photon energy Ex ~ 90 
eV. 

The only method of the attosecond pulse characteriza- 
tion, realized up to now, has been the so-called attosec- 
ond streak camera 0, 11111111, • The attosecond 
pulses under investigation ionize atoms of a gas target in 
the field of a co-propagating IR laser pulse. The bunch 
of photoelectrons presenting a replica of the XUV pulse 
interacts with the field of the few-cycle IR pulse. The 
streaking effect is achieved by measuring the photoelec- 
tron yield at different delays between the XUV and laser 
pulses. The resulting energy and angular distributions of 
photoelectrons, bearing information on the original at- 
tosecond XUV pulse, are measured by methods of elec- 
tron spectroscopy. Several theoretical methods are used 
to extract this information from the photoelectron spec- 
tra. They are based on the classical equation of motion 
0, the semiclassical approach HQ or quantum mechani- 



cal calculations within a strong field approximation |7j,|a|. 
In this way, it has been proved that single pulses with a 
duration of as low as 250 as H and even 150 as have 
been produced. The intensity of the XUV pulse has been 
also measured Q. In addition, the chirp was demon- 
strated as being close to zero [3j. In spite of the obvious 
success of the streak camera method, there are still many 
problems in understanding the detailed dynamics of the 
laser-dressed photoionization of atoms by the ultra-short 
pulses, the accuracy and limitations of the method. The 
latter problem has been recently discussed in detail in 
Ref. H3- 

In this Letter we present fully quantum mechanical cal- 
culations of the double differential cross section of atomic 
photoionization by an ultrashort XUV pulse in the pres- 
ence of a few-cycle intense IR pulse. The results provide 
complete information about the process which may be 
used either for comparison with an experiment or as a 
test-ground for approximate models. We present also a 
simple model which permits us to produce a quick analy- 
sis of the experimental data and to study the sensitivity 
of the measurements to the parameters of the attosecond 
pulses. 

The photoelectron spectra and angular distributions 
are calculated by solving the time-dependent Schrodinger 
equation numerically The atom is described by the 
independent-electron model, and we assume that the or- 
bitals of all electrons except the active one are frozen. 
This approximation is valid for not very strong laser fields 
( £l <C 1 a.u.). We suppose that both IR and XUV pho- 
tons are linearly polarized and their polarization vectors 
are parallel. In this case there is an axial symmetry with 
respect to the axis of polarizations which we choose as 
z-axis. By expanding the active electron wave function 
in spherical harmonics 
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i>(r,t)=J2u t (r,t)Y m (9, V ) (1) 
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and substituting the expansion into the time-dependent 
Schrodinger equation we obtain the following set of equa- 
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tions for the coefficient functions ut(r,t) 
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Here and in the following atomic units are used unless 
otherwise indicated. In Eq. (J2J) Ui(r) is a single elec- 
tron potential, C{£, £') is the angular part of the dipolc 
matrix element, £o is the binding energy of the active 
electron, is an electric field of the laser pulse while 

£x(t) is an envelope of the XUV pulse. Since the XUV 
field is weak, the first order perturbation theory and the 
rotating wave approximation are applied. Then in the 
dipole approximation, only transitions £q — > to ± 1 are 
possible. In contrast, the laser IR field is strong, it mixes 
the states with different orbital angular momenta. In or- 
der to achieve sufficient accuracy, the partial waves up to 
Lmax ^ 40 — 60 have been included. 

The system of equations has been solved numeri- 
cally using the split-propagation method. The details of 
the numerical procedure will be published elsewhere (see 
also Refs. julll3)- Calculating the set of the functions 
ug(r,t) at t — ► oo we have evaluated the partial ampli- 
tudes of photoionization by projecting the functions onto 
the corresponding continuum functions x~_ (E;r): 

/>oo 

Ae(E) = exp[irit(E)] / dru e {r,t +oc) X W (£; r), (3) 
Jo 

where pAE) is the photoionization phase (for details see 
Ref. 13]). Knowing the amplitudes, Ae(E), one can 
evaluate the double differential cross section (DDCS) 
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Here a is the fine-structure constant and K — 
J dt\£~x{t)\ 2 . By integrating Eq. Q over the emission 
angle one can obtain the photoelectron spectrum. 

As an example we consider the 3s subshell photoion- 
ization in Ar. The atom is described within the Hartree- 
Slater approximation with the same self-consistent 
potential U(r) for each £. The laser pulse is assumed to 
have the form 



£ L (t) = £ L (t)cos [uj L (t -t l ) + 4> l 
with the envelope given by the expression 
1 
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where tl gives the full width at half maximum (FWHM) 
and 4>l is the carrier-envelope phase. We set (j>L = 
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FIG. 1: (Color online) A 2D plot of the Ar(3s) photoelec- 
tron spectra (electron energy - x-axis) at 0° < i? < 15° as a 
function of the time delay in a.u. (1 a.u. = 24.2 as) between 
the laser and the XUV pulses (y-axis). The calculations were 
performed for the IR pulse of the cosine-type with the dura- 
tion t l = 5.0 fs, the laser intensity is 10 13 W/cm 2 . The XUV 
pulse duration is rx = 250 as and the carrier frequency is 90 
eV. The calculated binding energy of Ar(3s) is 28.6 eV 



(cosine- type pulse). The XUV pulse is delayed with re- 
spect to the laser pulse by td 



£ x [t) =£ x (t- t d ) cos [cox (t - t d )} 



(7) 



For the envelope £~x{t) we have chosen the hyperbolic 
secant shape. The parameters of the laser and the XUV 
pulses are shown in the caption to Fig. 1. This figure 
shows the two-dimensional (2D) plot of the calculated 
spectra as a function of the delay between the laser and 
the XUV pulses. The results shown in Fig. 1 are typical 
for the streaking experiment |2j, |3|, |4| . 

If the duration of the XUV pulse is much shorter than 
the period of the laser field, Tx < Ti, then one can con- 
sider the whole process as a two-step one |(| : first step is 
photoionization of the atom by the XUV pulse which is 
independent of the laser field, second step is a transport 
of the emitted electron by the laser field to the final state. 
Within the linear response approximation, the cross sec- 
tion of photoionization by a short pulse (without the IR 
field) can be presented as 01 



dp Q 



= ^== \D(E Q )\ 2 \f(E + \e \ uxt iftOM 



(8) 

Here po is the linear momentum of the photoelectron im- 
mediately after photoionization (# is the corresponding- 
emission polar angle, E is the electron energy), D(E ) 
is the dipole matrix element for transition to the con- 
tinuum state with the energy Eq, /(e) is the Fourier 
transform of the XUV pulse at the energy Eq + \eo\, and 
a = 2-KuoxotK~ 1 . The factor Y^tfo) describes the angu- 
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lar distribution of photoelectrons for the particular case 
of the s— >p transition considered here. 

With the IR field, the cross section modifies as follows 
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where p is the final momentum of the electron, 
is the transformation determinant. Now we are to link 
the variables ($oi Eq) with the experimentally observable 
quantities E). The basic relations for the components 
of the momentum p = {p\\,p±} parallel and perpendicu- 
lar to the field are 



P\\ =Po\\ 



from which it follows 



£ L (t) dt= Po{] -A l (t), (10) 
p±=Po±, (11) 



d(po) 



d(p) 



1 and 



sin $ = V ' E/Eq sin § . 



(12) 



Here r is the time of the photoelectron emission and 
Al(t) is a parallel component of the vector potential of 
the laser field in the Coulomb gauge. From Eqs. i|l(J|) 
and l|12fl it follows that 

E a = E + V2EA l (t)cos$ + A\{t)/2. (13) 

The final approximate expression for the DDCS reads: 
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where E ,i) and E,i!} are connected by Eqs. l(T^l[T^|) . 

In Fig. 2(a) we compare the DDCS calculated by the 
numerical solution of the Schrodinger equation and by 
using Eq. (|14|> in the case when the delay time is 7 fs. 
The agreement between the two calculations is almost 
perfect. Note that at the chosen delay ft^d) is almost 
zero. Similar agreement is obtained for all other delays 
when the laser field is close to zero. In contrast, for the 
case oi td — 6.5 fs shown in Fig. 2(b), the distribution 
given by the exact computation is substantially broader 
than that given by the approximate formula. Similar 
discrepancy is obtained for other delays when the value 
of the laser field is close to its extrema. 

From the angular distributions calculated for each time 
delay we obtain the photoelectron energy spectra at a 
given angle by integrating over some acceptance angle 
(15°) of the detector. These spectra are fitted with the 
Gaussian function A exp[— [E — E) 2 /T 2 ] in order to de- 
termine their widths. Variation of the spectral width 
with the time delay is commonly used in order to deter- 
mine the parameters of the XUV pulses 0,0. In Fig. 
3 we show the widths T of the calculated spectra as a 
function of the delay time for two angles $ — 0° and 
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FIG. 2: (Color online) 2D plots of the calculated DDCS for 
Ar(3s) photoionization for different time delays between the 
laser and the XUV pulses: (a) t d = 290 a.u. (7 fs), (b) 
td — 270 a.u. (6.5 fs). The parameters of the pulses are the 
same as in Fig. 1. Contour plots show the results obtained 
by solving the Schrodinger equation; color (grey) scale plots 
show the results of model calculations using Eq. fl2). The 
minimal and maximal values are the same for both plots. 



180°. Comparing the results of the exact calculations 
(symbols) with the calculations using the approximate 
formula Eq. (|14|) (dotted curves) we have noticed that 
the difference is always proportional to the value of the 
IR laser electric field at the moment of the photoelectron 
emission. As is known Q , the broadening of the electron 
spectra (for electrons emitted in the same direction) is 
determined by two factors. First, the spectrum becomes 
broader since the emitted electrons have different ener- 
gies due to the large spectral width of the XUV pulse. 
This type of broadening is taken into account by Eqs. 
(|13fl and (|14|) . Second, the electrons are emitted at dif- 
ferent times within the XUV pulse duration and, there- 
fore, at slightly different phases of the IR field that leads 
to an additional spread of the electron energies. Since 
tx -C Tl, this spread depends linearly on the field at 
the moment of electron emission, and in the first approx- 
imation it is proportional to the duration of the XUV 
pulse. Thus we suggest a semiempirical expression for 
the width: 



F — F 



/ 2^f L (^)r 1/2 /2 : 



(15) 



where r m is a model width, calculated from the spec- 
tra obtained using Eq. (fHfjl and t\ii is a FWHM of the 
XUV field envelope. Using this expression we obtained 
the values, shown by solid lines in Fig. 3, which are in 
excellent agreement with the results of the exact calcula- 
tions shown by open symbols, circles (forward emission) 
and squares (backward emission). 

Besides very good agreement between the calculations, 
an interesting feature can be seen in Fig. 3, namely a 
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FIG. 3: (Color online) The width of the photoelectron spec- 
trum as a function of the delay time for two XUV pulse du- 
rations: (a) tx = 100 as, (b) tx = 250 as. Black sym- 
bols and curves correspond to forward emission # = 0°, blue 
(gray) symbols and curves correspond to backward emission 
■d = 180°. Exact calculations: circles (d — 0°) and squares 
($ = 180°). Dotted curves show Y m obtained by model cal- 
culations. Solid curves show the results obtained by Eq. 1151 . 
All parameters in the calculations are the same as in Fig. 1. 



double- maximum structure in the resulting model curve 
for the pulse duration of 100 as. This is not an artefact 
since the exact calculations with small steps in time de- 
lay around td — 300 a.u. perfectly agree with the model 
calculations. The structure is due to the fact that at the 
zero-crossing points where £r,(~td) = 0, the width should 
be equal to that obtained from Eq. 1)140. and therefore 
a minimum is formed. The depth of the minimum is de- 
termined by the contribution of the pulse spectral width 
(dotted curve). For longer pulses (250 as) this contribu- 
tion is smaller and the minimum is deeper. An analysis 
of the results obtained with model expression Eq. 114|) 
shows that for very short XUV pulses ( < 100 as) the 
model describes the width quite well, as well as describ- 
ing the whole structure of the DDCSs. Thus, for such 
short pulses one can keep track only on the spectrum of 
the XUV pulse but not on its duration and the chirp. 
With an increase of the XUV pulse duration, the first 
term in Eq. (|15fl is decreasing while the second one is in- 
creasing and becomes dominant. The second term gives 
information on the FWHM of the XUV envelope. The 
latter quantity can also depend on the chirp, although 
this dependence is not strong. Thus, the chirp can be 
obtained only for the pulses for which the second term 
in Eq. I|15|) is sufficiently large. The value of this term 



is proportional to the strength of the IR field and to the 
duration of the XUV pulse. Thus stronger IR fields and 
longer XUV pulse durations are favorable for chirp reg- 
istration. 

In conclusion, we have developed a fully quantum me- 
chanical theory of the streaking measurements in the at- 
tosecond region based on numerical solving of the time- 
dependent Schrodinger equation which includes the elec- 
tron interaction with the atomic core as well as with 
the XUV and the laser fields. We suggest also a sim- 
ple model expression for the width of the photoelectron 
spectra which agrees very well with the exact calcula- 
tions. Using this expression it is easy to analyze the 
sensitivity of the width to the parameters of the XUV 
pulse. 
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